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Abstract. The study of thermal emission from isolated, cooling neutron stars plays a key role in probing the 
physical conditions of both the star crust and the core. The comparison of theoretical models for the star thermal 
evolution with the surface temperature derived from X-ray observations of sources of different age is one of the 
main tools to investigate the properties of the interior and constrain the equation of state. Here we propose 
to use population synthesis studies as an independent approach to test the physics governing the star cooling. 
Theoretical Log N-Log S distributions depend on the assumed neutron star thermal evolution. We have computed 
distributions for several different cooling scenarios and found that comparison with the observed Log N-Log S 
of isolated neutron stars is effective in discriminating among cooling models. Among the eleven cooling models 
considered in this paper, all of which may reproduce the observed temperature vs. age diagram, only at most 
three can explain the Log N-Log S distribution of close-by cooling neutron stars. The Log N-Log S test, being a 
"global" one and despite some limitations, appears indeed capable of complementing the standard temperature 
vs. age test used up to now. 
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1. Introduction 

The determination of the equation of state (EOS) at or 
above nuclear densities is a long sought goal in high en- 
ergy astrophysics. In this respect direct observations of 
neutron stars (NSs) may provide invaluable insight into 
several key issues of fundamental physics, like quantum 
chromodynamics/electrodynamics, superfluidity and su- 
perconductivity, that could not be otherwise tested under 
laboratory conditions. An ideal way to place tight con- 
straints on the EOS is the simulta neous mass and ra- 
dius measurement of NSs (see e.g. iLattimer fc PrakashI 
l200l|) . Precise mass determinations have been obtained 
for NSs in binaries, especially in radio pulsar systems (e.g. 
iThorsett fc Chakrabatvll999tl . Simultaneous mass and ra- 
dius measurements are presently available for a few X-ray 
binaries, although masses derived from gravitational red- 
shift of spectral line s are still uncertain in many cases (e.g. 
ICottam et al]l2002|) . Mass estimates can be obtained also 
from quasi-periodic oscillation measurements ijBulik et al.l 
but results based on this, as well as other ap- 
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proaches, are still quite model dependent. The study of 
glitches observed in radio pulsars, and rec ently in anoma- 
lous X-ray pulsars (A XPs) as well (e.g. iHorvathI l2004t 
iDairOsso et"al]l200.'j) . promise further insight into the 
understanding of the internal structure of NSs, as well 
as future observations of neutrino and gravitational wave 
emission from neutron star sources. 

For the time being, however. X-ray observations of NSs 
are, and will still be for some time, central for our under- 
standing of the star interior. Isolated neutron stars that 
emit at X-ray energies as they cool are particularly promis- 
ing in this respect. Their thermal radiation directly comes 
from the star surface, carrying information on the physical 
conditions of the emitting matter, in particular on the star 
surface temperature. Thermal X-ray emission has been de- 
tected from about 20 isolated NSs so far, including normal 
radio pulsars, central compact objects in supernova rem- 
nants (CCOs in SN Rs), radioquiet NSs, AXPs and soft 
7-repcatcrs (sec e.g. IPavlov fc ZavlinI l2003t iKaspi et alJ 
12004: Haberl 2004 for reviews). 

Significant progress in the understanding of NS ther- 
mal evolution has been made in recent years and 
cooling curves have been computed by several groups 
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(see e.g. jK aminker ct al.' '2002'; 'Ts uruta etal] l2002t 
IPaee et aP 12004: Blaschke ct al. 20041 and references 
therein). The present work is based on NS cooling calcu- 
lations performed in the Nuclear medium cooling scenario 
l|Blaschkc ct al. 2004) which differs from the other above- 
mentioned approaches in a consistent inclusion of medium 
effects. Processes of internal heating l|Tsuruta et al.lE'002l) 
are not included. Since the cooling history crucially de- 
pends on the assumed physical conditions inside the star, 
comparison with observations may rule out some models 
in favor of others. To study this we will vary assump- 
tions on the nuclear pairing gaps, the relation between 
crust and surface temperature as well as presence or ab- 
sence of pion condensation. A customary way of testing 
predictions of cooling calculations is to construct a tem- 
perature vs. age (T-t for short) plot for the largest sample 
of sources. Despite its wide application and undisputed 
usefulness, this test has a number of limitations. In this 
paper we suggest use of the LogN-LogS distribution of 
close-by NSs as an additional probe for NS cooling mod- 
els. The idea is based on the comparison of present obser- 
vational data with NS population synthesis calculations 
in which cooling curves are one of the ingredients. Our 
approach extends the actual calculation of observational 
properties of the population o f close-by yo ung NSs that 
has already been developed bv lPopov et al., (2003, 2005) 
(hereafter Paper I and II) by including nuclear medium 
effects in the cooling code and by investigating the con- 
tribution of massive progenitors in the Gould Belt to the 
local NS population in more details. 

The paper is organized as follows. In section El we dis- 
cuss the advantages and limitations of the two methods in 
extracting information from observational data. Section |31 
presents our population synthesis model in more detail, es- 
pecially concerning the choice of the set of cooling curves 
which will be used in our calculation. We present our re- 
sults and discuss them in sectional Scction[3contains our 
conclusions. 

2. The two tests 

In this section we briefly discuss and compare the ca- 
pabilities of the conventional T-t test and of the LogN- 
Log S test we propose here. Our main conclusion is that 
the two approaches should be used together as their ad- 
vantages and disadvantages are mostly complementary to 
each other especially in the case of available samples of 
observed objects. 

2.1. The T-t test 

The T-t test is the most appropriate one to compare re- 
sults of thermal evolution calculations with observations. 
An important advantage of this test is that there are 
no additional theoretical uncertainties except those con- 
nected with the cooling model: theoretical cooling curves 
do not depend on unknown (or poorly known) astrophys- 
ical parameters but only on the input physics of the star 



interior (for a discussion of co oling processes and r efer- 
ences to earlier pape rs, see e.g. lYakovlev et al .1179991 and 
iBlaschke et"al] l2004j) . Still, there are some well known 
drawbacks to this method. 

For sources associated with radiopulsars or CCOs, the 
star age is usually estimated from the spindown age P/ P, 
or inferred from the age of the supernova remnant. To 
which extent these determinations are indeed representa- 
tive of the neutron star age is still uncertain. The situa- 
tion is even worse for ROSAT isolated neutron stars which 
are not associated with a SNR and are with no exception 
radio-silent. Current age estimates for the two brightest 
objects in this class (RX J1856. 5-3754 and RX J0720.4- 
3125), based on dynamical considerations, should be re- 
garded only as guesses. 

The T-t test is not very sensitive to objects with ages 
^ 10^ yrs. There are two reasons for this. First, there 
are only a few sources older than this value to which the 
test can be applied. Then, cooling curves sharply drop at 
the photon cooling stage. Shifts between different cooling 
curves are comparable with data error bars. So, it is diffi- 
cult to discriminate between models that differ mainly in 
this respect. 

Cooling calculations provide the star temperature at 
the core-crust boundary and the actual surface temper- 
ature is then obtained applying a bridging formula (e.g. 
lTsurutalll979t lYakovlev et al.ll2004|) . Detailed modehng of 
heat transport in the highly magnetized envelope indicates 
that the surface temperature may be influenced by several 
effects, among which the magnetic fie ld distribution inside 
the star is very important ijCeppert et al.-2004 'l . The tem- 
perature is derived from a spectral fit to the data and, as 
such, depends on the assumed emission model for the sur- 
face. Different models (blackbody, H/heavy elements at- 
mosphere with/without magnetic field, solid surface) give 
values which can differ by a factor of a few. 

Recently, it has been suggested that the Temperature- 
Age test should be refined by the introduction of a bright- 
ness constraint (BC), which demands that cooling curves 
for objects within a mass interval anticipated as typical 
should not appear significantly brighter than t he bright- 
est of already observed objects for a given age ijCrigoriaiJ 
I2OO5I hereafter G05). It has been shown in this work that 
some of the cooling models discussed in BGV which would 
violate the BC become acceptable when the crustal prop- 
erties are modified. These modifications do not affect old 
stars with ages 10^ yrs for which the LogN-LogS test 
is relevant. Due to our lack of knowledge of the mass dis- 
tribution and atmospheric properties of cooling compact 
stars, we must be cautious in applying the strong BC. 
Therefore we refer in this paper to results of the T-t test 
with and without BC. 

Finally, the data sample is non-uniform: there are 
sources of different types, and, more important, the sam- 
ple is neither flux, nor volume limited, and is a strongly 
selected one, containing objects for which age and tem- 
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perature estimates are available^. Clearly these sources 
are not representative of any real population of NSs. 

2.2. 7/76 Log N-Log S test 

The Log N-Log S distribution is a widely used tool in many 
branches of astronomy. For isolated NSs such an approach 
has been already used by iNeuhauser h Triimperl l)l999|) 
and IPodov et al.l l|2000l) to probe the origin of isolated 
NSs in the solar proximity. The main conclusion of these 
investigations is that the observed Log N-Log S cannot be 
easily explained assuming that the local population of NSs 
originated only in the Galactic disc. As shown in Paper 
I and n, accounting for massive progenitors in the Gould 
Belt reconciles theoretical predictions with the data. 

One immediate advantage of the Log N-Log S test is 
that, at variance with the T-t test, no degree of arbitrari- 
ness is introduced when observational data are analyzed: 
both the fluxes and (of course) the number of sources are 
well measured. In addition, this approach is a "global" 
one. In our scenario it would not be possible to explain 
some particular sources by invoking slight changes in the 
cooling physics. Once the parameters of the model other 
than those related to the cooling process are fixed (see 
section a particular cooling curve either fits the pop- 
ulation as a whole or not. Furthermore, the Log N-Log S 
sample is a uniform one, i.e. objects are flux (and prob- 
ably volume) limited, and no strong selection criteria are 
introduced. 

For the Log N-Log S test the only necessary observa- 
tional piece of information is the ROSAT count rate. The 
method can be applied to objects with unknown ages. This 
makes it possible to include, for example, all the ROSAT 
X-ray dim NSs, and 3EG J1835+5918 (the Geminga twin) 
in the testing sample. The Log N-Log S test is mostly sen- 
sitive to NSs older than ~ 10^ yrs. Older sources dominate 
in number, and in the solar proximity there are about a 
dozen of them in comparison to very few with t ^ 10^ 
years. ^ 

Nevertheless there are significant limitations too. One 
source of uncertainty is our incomplete knowledge of some 
important ingredients of the population synthesis model. 
These are discussed in more detail in the next section and 
concern the spatial distribution of the NS progenitors, the 
NS mass and velocity spectrum, and their emission prop- 
erties. However, all these issues, to some extent, can be 

^ Note that there are some young cooling NSs for which their 
age is not known, in particular those not observed as active 
radiopulsars. These sources can not be included in the sample 
of the Temperature- Age test. 

^ Of course, the Log N-Log S test can be sensitive to young 
objects if another uniform sample is used. For example, it is 
very important to make a population synthesis of sources in 
SN remnants. Then, as we mentioned above, differences be- 
tween cooling curves for large ages are tiny. However, in the 
Log N-Log S approach we, in some sense, integrate these small 
differences along the curves, so their impact on the final Log N- 
LogS distribution is significant. 



addressed by considering different cases believed to cover 
the entire range of acceptable scenarios. In addition, there 
could be unknown correlations among some of the quanti- 
ties we use to parametrize our model, so that they should 
not be treated as independent ones. Examples of such pos- 
sible correlations are those between the star kick and the 
internal structure (beca use of quark deconfinement, see 
iBombaci fc Popovll2004(l . and between the s tar mass and 
the magnetic field (beca use of fallback, see IPopov et alJ 
l2002t iHeger et al.l 120051) . A more severe problem arises 
in connection with the low statistics of the sample, since 
there are only about 20 thermally emitting NSs known to 
date. This implies that the bright end of the Log N-Log S 
relation comprises very few objects so that it is difficult to 
account for statistical fluctuations. We do not know much 
either about the properties of very faint sources, i.e. the 
dim end of the Log N-Log S distribution. 

3. The population synthesis model 

The main physical ingredients that enter our population 
synthesis model are: 

— the initial NS spatial distribution; 

— the kick velocity distribution; 

— the NS mass spectrum; 

— the cooling curves; 

— the surface emission; 

— the interstellar absorption. 

The calculation of the NS spatial evolution as they 
move in the Galactic gravitational potential follows that 
presented in Papers I and II . The same treatment of the 
interstellar absorption is retained and the kick distribution 
is that proposed by Arzoumanian et al. (2002). We do not 
account for atmospheric reprocessing of thermal radiation, 
and assume that the emitted spectrum is a pure black- 
body. Although this is clearly an oversimplification, it is 
a reasonable starting assumption and will serve for our, 
mainly illustrative, purposes. A more detailed description 
of surface emission may be easily accommodated in our 
model later on. For the time being, we perform our cal- 
culations for nine different sets of cooli ng curves among 
those discussed in Blaschk e et al ] <l2004 . hereafter BGV, 
to which we refer for all details. This issue, together with 
the initial spatial distribution of NSs and their mass spec- 
trum, is further discussed below. 

3.1. The initial NS spatial distribution 

Following the results of previous investigations (paper I), 
we take as an established fact that the population of 
nearby NSs is genetically related to the Gould Belt. The 
contribution of the Belt has dominated the production of 
compact remnan ts in the sola r proximity over the past 
~ 30 Myrs (see |Poppe]|ll997l for a detailed description 
of the Belt structure). About two thirds of massive stars 
in the ~ 600 pc around the Sun belong to the Belt 
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ifTorra et al.ll2000|) . Since they are bright objects, it is 
possible to track their positions exactly for example by 
using HIPPARCOS data. However, here we use a simpli- 
fied distribution of the progenitor stars, as discussed in 
Papers I and II. NSs are born in the Belt, for which a sim- 
ple ring-shaped structure is assumed, and in the Galactic 
disc. The outer belt radius is a model parameter and two 
values, 300 pc and 500 pc, were used, the second being an 
upper limit. The supernova rate was taken from ICrenieJ 
1(200 0^. and appears in agreernent w ith the historical rate 
estimated bv iTammann et al 

3.2. The NS mass spectrum 

Since cooling curves are strongly dependent on the star 
mass, the mass spectrum is one of the most important in- 
gredients and, unfortunately, one of the lesser known. We 
cannot r ely on the mass measurements in binary radio pul- 
sars fe.g. lThorsett fc Chakrabatvll999tl because they refer 
to a "twice selected" population (i.e. selection effects due 
to evolution in a binary can be important together with 
possible conditions necessary for radio pulsar formation). 
Probably not all NSs go through the active radio pulsar 
stage fe.g. lGotthelf fc Vasis ht 2000). and the properties of 
NSs in bi naries may be different to those of isolated ob- 
jects fe.g. IPodsiadlow ski et al.l l200,'^ . Note that the local 
NS mass spectrum can be different to the global NS mass 
spectrum in the Galaxy. Even stronger deviations can be 
expected between the mass spectra of local NSs and of 
those sources usually used for the T-t plot. 

As the population of NSs in 1 kpc around the Sun 
may be slightly different from the average galactic pop- 
ulation, we estimated the mass spectrum for these ob- 
jects directly (see Paper II for more details). The ba- 
sic idea is to use HIPPARCOS data on massive stars 
around the S un in conjunction with the calculations by 
IWooslev et al . (2002). Knowing the mass distribution of 
progenitors through their spectra l classes, we use a fit to 
a plot from lWooslev et'al] l|2002|) in order to obtain the 
NS mass from the mass of the progenitor. 

We use eight mass bins centered at M/Mq = 1.1, 1.25, 
1.32, 1.4, 1.48, 1.6, 1.7, 1.76. The adopted mass spec- 
trum is shown in Fig. ^ The lower limit for th e NS mass 
is still an open question. iTimmes et al.l ljl996|) suggested 
that there are no NSs with M <, 1. 27 Mrr,, although their 
conclusion is not definite (see also IWoosle v et al. 2002). 
For this reason we decide to use also a truncated mass 
spectrum, in which the first bin is suppressed and all ob- 
jects originally contained there are added to the second 
one. Each bin corresponds to one of the calculated cool- 
ing curves. According to the mass spectrum, each curve 
has a statistical weight of 31.75%, 25.75%, 11%, 28.125%, 
0.875%, 1.125%, 0.75%, and 0.625%. For the truncated 
one the weights of the first two bins are replaced by 0.0 
and 57.5% respectively. We note that sampling also rela- 
tively low masses is important since low-mass NSs seem 
to be required to interpret data on the T-t plot (BGV). 




2.0 



Fig. 1. The adopted mass spectrum, binned over eight 
intervals of different widths. 



3.3. Cooling curves 

In their recent paper iBlaschke et al presented six- 

teen sets of cooling curves. Each set contains models for 
several values of the star mass while different sets refer 
to different assumptions on heat transport in the crust 
and on the physical processes in the NS interior. Five 
of these models are unable to reproduce the observed 
temperature-age plot and will not be considered further. 
From the remaining eleven sets, all of which give results 
not in contradiction with observations (see, however, the 
discussion in sec.0J, we select nine representative models 
for our population synthesis calculations (models I - IX 
in Table We add two models (X and XI) from the re- 
cent analysis in G05, which correspond to models HI and 
VII, respectively, when calculated with different crustal 
properties. All of them have superfluid nuclear matter and 
medium modifications of the neutrino processes. They dif- 
fer in the assumptions about the superfluid gaps, the pres- 
ence/absence of a pion condensate and the properties of 
the neutron star crust. The latter governs the relationship 
between the temperature of outermost core layer (Tin) to 
that of the star surface (Tg). The main characteristics of 
the selected models are summarized in Table Q All of the 
eleven models satisfy the temperature - age test according 
to BGV, whereas only six of them fulfill the additional 
brightness constraint introduced in G05. The last column 
anticipates the results discussed in the following section 
and shows if the model complies with the LogN-LogS 
test. 



4. Results and discussion 

In each run we calculate 5000 individual tracks for the spa- 
tial evolution of a single star with a time step of 10^ yrs. 
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Fig. 2. Model I. Left: cooling curves for (from top to bottom) 1.1, 1.25, 1.32, 1.4, 1.48, 1.6, 1.7, 1.76 M©. Right: the 
corresponding Log N-Log S distribution for Rbcit — 300 pc and non-truncated mass spectrum (full line) and 500 pc 
and truncated mass spectrum (dotted line). See text for details. 




Log t (years) Log S (cts/s) 

Fig. 3. Same as in Fig. |21 for Model IIL The dashed line refers to a calculation in which the full (non-truncated) 
mass spectrum was used and Rbdt was assumed to be 500 pc. Model X should produce nearly the same LogN-LogS 
distribution as it differs only in the type of the crust. 



Each track is applied to all eight (or seven for the trun- 
cated spectrum) masses, and the thermal evolution is fol- 
lowed by the correspondig cooling curve. Results are then 
collected according to the statistical weight of each mass 
bin. 

Results are summarized in Tableland Figs. l2l7l which 
refer to a selected sub-sample of the cooling curve sets 
listed in that Table. In the left panels of Figs. I2I7I the 
corresponding cooling curves for the various masses are 



shown. Results are plotted for ages > 10** yrs and temper- 
atures above 10^ K. The right panels illustrate the LogN- 
Log S distributions computed for the same sets of cooling 
curves. All models have been calculated using both the full 
mass spectrum and the truncated one, although results for 
the latter are not shown in all cases (see the following dis- 
cussion). Two values of the outer radius of the Gould Belt 
have been used to test the dependence of our calculation 
on the assumed geometry, i?boit — 300 pc and 500 pc. 
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Log t (years) Log S (cts/s) 

Fig. 4. Same as in Fig.Ofor Model VI. Model IV gives quite similar results since the change of the crust model from 
C to E does not affect the log N - log S distribution. 




Log t (years) Log S (cts/s) 

Fig. 5. Same as in Fig. O for Model VII. Model XI should produce nearly the same LogN-LogS distribution as it 
differs only in the type of the crust. 



Theoretical distributions are superimposed on the ob- 
served LogN-LogS for isolated NSs. The data points 
are derived from the sample of thirteen sources listed 
in Paper I. Error bars correspond to Poisson statis- 
tics and are plotted to illustrate the statistical signifi- 
cance of the points. We note, however, that there can 
be more unidentified sources, especially at fluxes below 
0.1 cts s~^ (see^utledgc et al. 2003 for a recent discus- 
sion). In this respect, the fact that the two points at 
the lowest fluxes lie below the general trend of the ob- 
served LogN-LogS is not surprising. In addition, the ob- 



servational upper limit derived bv lR,utledge et all l(200.'^ 
on the number of fainter sources is also shown (marked 
as BSC in the figures). Bounds on the total number of 
s ources with flux > .2 cts s~^ have been also presented by 
iSchwope et all l|l999>) . on the basis of the ROSAT Bright 
Sources Catalogue (BSC). 

The comparison of the predicted and observed Log N- 
LogS distributions in Figs. [2171 indicates that at most 
three cooling models (model I and possibly models VIII 
and IX) are in agreement with the data. All the others 
substantially overpredict the observed number of sources 
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Fig. 6. Same as in Fig. El for Model VIII. 
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Fig. 7. Same as in Fig. El for Model IX. 



at all fluxes, even though they comply with the T-t test. 
This latter statement deserves further comment. BGV did 
not reject models II- VII on the basis of the fact that the 
cooling curves cover, for the assumed mass range, the en- 
tire region in the T-t plane where the observed sources lie. 
This approach is sound, and it represents the only possible 
option to discriminate among different cooling scenarios at 
the zero level, i.e without introducing additional informa- 
tion. Clearly, if some assumptions on the NS mass distri- 
bution are made, the T-t test can be used to exclude some 
further models. If the same mass distribution discussed in 
Sec. 13.21 is applied to BGV sets of cooling curves, one is 
immediately led to discard models II-V and VII because 



they predict quite high temperatures for low mass NSs 
(M ^ I.3M0) which, according to our mass spectrum, 
are very abundant. Although such hot objects would be 
detectable even to large distances they are not actually 
observed (this issue is further discussed in G05, see also 
below). However, care must be taken in using such an ar- 
gument. Our mass spectrum is meant to be representative 
of the local population of isolated NSs, and its application 
to the very biased and limited sample of objects which 
can be placed on the T-t plots is uncertain. Nevertheless, 
one can reverse the argument by saying that the two tests, 
when provided with the same amount of information, yield 
results that are broadly consistent, as they should. Still, 
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we note that a strict interpretation of the LogN- LogS 
test results in the exclusion of two further models (VI and 
VIII) and that even in the most optimistic case model VI 
is rejected. 

The Log N-Log S test is not equally sensitive to changes 
in the three main groups of parameters in Table ^ (pres- 
ence or absence of the pion condensate, gaps, and type of 
the crust). Mostly the test reacts to changes in the gap 
parameters. Conversely, changes in the type of the crust 
are not very important as here we consider a sample of rel- 
atively old NSs (nevertheless uniform samples of younger 
NSs can also be studied). Variations in the crust proper- 
ties are discussed in detail in G05, where the brightness 
constraint was introduced. G05 shows that models II-V 
and VII can be ruled out when the absence of very bright 
young NSs is considered as a constraint. The fact that the 
Log N-Log S test leads to the same conclusion is however a 
completely independent result, since objects from different 
age ranges are used in the two approaches. 

As mentioned earlier, the Log N-Log S test applied to 
close-by NSs is not sensitive to changes in the crustal prop- 
erties. This implies that if a cooling model with a given 
crust is rejected then its variants with other types of crust 
(at least from the set considered here) can be ruled out, 
too. For example, additional models X and XI (see Table 
^ do not satisfy the test because their twins (models III 
and VII, respectively) do not (see Figs. El and 01. In G05 
it has been shown that it is possible to fulfil the BC by 
changing the crust. For example, model III with crust E 
can fit the data. However, on the basis of the Log N-Log S 
test this model is ruled out. This is because the latter 
test mainly "feels" the changes in gaps. And if they are 
"wrong" then changes in the crust properties are of no 
avail. 

Jumping to the conclusion that the present analysis 
can provide direct information on the physical state of 
star interior would be premature. Our adopted scenario 
contains a number of uncertainties, as discussed in sec- 
tions |21 and |2| Nevertheless, we believe that the case pre- 
sented here convincingly shows that by combining theoret- 
ical cooling curves with populations synthesis calculations 
one has the potential to discriminate between competing 
cooling scenarios. The main outcome of this investigation 
is not that models I, and possibly VIII and IX, fit the data 
while others do not. This may be the result of our start- 
ing assumptions. What matters is the fact that, within the 
same set of assumptions, different cooling models produce 
different results when compared with observations. 

Current limitations of this approach are due to our 
present incomplete knowledge of some key issues, chiefly 
the NS mass spectrum, their surface emission properties 
and initial spatial distribution. We attempted to account 
for some of these uncertainties in our model by consider- 
ing different configurations which should bracket the true 
behavior. We note that in all the six cases for which no 
agreement has been found between theoretical and ob- 
served distributions, this is largely independent of the as- 
sumptions we introduced. Using different values of i?boit 



or taking variants of the mass spectrum does not help in 
reconciling predictions with the data. However, we caution 
that other effects, like those introduced by the proper in- 
clusion of an atmospheric model, may be important. This, 
and other issues, will be the subject of future work. 

In all the cases we examined (see again Figs. l2l7|) . with 
the possible exception of model VIII, the capability of 
our test to discriminate between different cooling scenar- 
ios seems to be quite robust. The three models (I, VIII 
and IX) which can reproduce the observed Log N-Log S 
for our choice of parameters have also been considered 
by BGV as the most realistic ones. Among models I-IX 
which according to BGV are not in contradiction with the 
Temperature-age test, the three mentioned above are the 
theoretically most appealing ones since either the super- 
fluid gaps were calculated with the same nucleon-nucleon 
interaction that formed the basis for the equation of state 
and thus the structure of the neutron star configurations 
(models I and IX) or the gaps were modelled such as to 
mimic these results (model VIII). 

5. Conclusions 

In this paper we suggest adopting the Log N-Log S test 
as an useful addition to the standard T-t test in probing 
neutron star cooling models. To illustrate the capabilities 
of the proposed ap proach, we app l ied it to nine sets of 
cooling curves from Blasch ke et al ](^ooi). Out of sixteen 
sets described in that p aper and two additional ones taken 
from ICrigorianl l)2005l) . these eleven produce results that 
are not in immediate contradiction with the T-t test. The 
application of the Log N-Log S test rules out at least eight 
out of eleven investigated cooling models, resulting in just 
three models able to pass both tests. Requiring that the 
tested cooling models should fulfill in addition the more 
stringent bightness constraint of G05, there are still six 
models left out of which only 50% pass the Log N-Log S 
test. 

One of the most challenging questions for the applica- 
tion of the test suggested in this paper is to use it next for a 
possible discrimination between purely hadronic compact 
star cooling scenarios and hybrid ones for stars haying a 
color superconducting quark core iGrigorian et~a 
which have already successfully passed the T-t test. Our 
conclusion is that the Log N-Log S may therefore become 
a powerful strategy in uncovering the properties of dense 
nuclear matter under the extreme conditions in neutron 
star interiors. 

Acknowledgements. We thank D. Yakovlev and other mem- 
bers of the St-Petersburg group for many helpful discussions. 
We are indebted to D. Voskresensky and M. Prokhorov for 
their constant interest and support. Comments by the ref- 
eree were very useful and helped to improve the paper. S.P. 
acknowleges a postdoctoral fellowship from the University of 
Padova where most of this work was carried out and a fellow- 
ship from the "Dynasty" Foundation. The work of S.P. was 
partly supported through RFBR grant 03-02-16068. H.G. is 
grateful to the Department of Physics, University of Padova, 



S. Popov et al.: Population synthesis and neutron star thermal evolution 



9 



for hospitality and acknowledges financial support from the 
Virtual Institute "Dense hadronic matter and QCD phase 
transition" of the Helmholtz Association (grant No. VH-VI- 
049) and from Deutsche Forschungsgemeinschaft (grant No. 
436 ARM 17/4/05). The work was partially supported by the 
Italian Ministry for Education, University and Research under 
grant PRIN-2002-027245. 



References 

Arzoumanian, Z., Chcrnoff, D.F., Cordes, J.M. 2002, ApJ 
568, 289 

Blaschke, D., Grigorian, H., Voskresensky, D. 2004, A&A, 

424, 979 (BGV) 
Bombaci, I., Popov, S.B. 2004, A&A, 424, 627 
Bulik, T., Kluzniak, W., Zhang, W. 2000, A&A, 361, 153 
Cottam, J., Paerels, F., Mendez, M. 2002, Nature, 420, 51 
DairOsso, S., Israel, G.L., Stella, L. et al. 2003, ApJ, 599, 

485 

Gotthelf, E.V., Vasisht, G. 2000, in Proc. lAU Coll. 177 
"Pulsar astronomy - 2000 and beyond" , eds. M. Kramer, 
N. Wex, & N. Wielebinski, ASP Conf. Series 202, 699 

Geppert, U., Kiiker, M., Page, D. 2004, A&A, 426, 267 

Grenier, LA. 2000, A&A, 364, L93 

Grigorian, H., Blaschke, D., Voskresensky, D.N. 2005, 

Phys. Rev. C, 71, 045801 
Grigorian, H. 2005, astro-ph/0507052 (G05) 
Haberl, F. 2004, Adv. Sp. Research, 33, 638 
Heger, A., Woosley, S.E., Spruit, H. 2005, ApJ, 626, 350 
Horvath, J.E. 2004, Int. J. Mod. Phys., D13, 1327 
Kaminker, A.D., Yakovlev, D.G., Gnedin, O.Yu. 2002, 

A&A, 383, 1076 
Kaspi, V.M., Roberts, M.S.E., Harding, A.K. 2004, to 

appear in "Compact Stellar X-ray Sources", W.H.G. 

Lewin and M. van der Klis Eds. |astro-ph/0402136| 
Lattimer, J.M., Prakash, M. 2001, ApJ, 550, 426 
Neuhauser, R., Triimper, J. 1999, A&A, 343, 151 
Page, D., Lattimer, J.M., Prakash, M. et al. 2004, ApJS, 

155, 623 

Pavlov, G.G., Zavlin, V.E. 2003, in proceedings of the XXI 
Symposium on Relativistic Astrophysics, R. Bandiera, 
R. Maiolino, F. Mannucci Eds. World Scientific 
PubHshing (Singapore), p. 319 astro-ph/0305435 

Podsiadlowski, P., Langer, N., Poelarends, A.J.T. et al. 
2004, ApJ, 612, 1044 

Poppel, W. 1997, Fund. Cosm. Phys., 18, 1 

Popov, S.B., Colpi, M., Prokhorov, M.E. et al. 2000, ApJ, 
544, L53 

Popov, S.B., et al. 2002, in Proc. of the Third International 
Sakharov Conference on Physics, ed. A. Semikhatov et 
al. , Scientific World: Moscow, p 420 astro-ph/0210688 

Popov, S.B., Colpi, M., Prokhorov, M.E. et al. 2003, A&A^ 
406, 111 (Paper I) 

Popov, S.B., TuroUa, R., Prokhorov, M.E. et al. 2005, 
Ap&SS, 299, 117 astro-ph/0305599 (Paper II) 

Rutledge, R.E., Fox, D.W., Bogosavljevic, M. et al. 2003, 
ApJ, 598, 458 



Schwope, A. D., Hasinger, G., Schwarz, R. et al. 1999, 

A&A, 341, L51 
Takatsuka, T., Tamagaki, R. 2004, Prog. Theor. Phys., 

112, 37 

Tammann, G.A., Loffler, W., Schroder, A. 1994, ApJS, 92, 
487 

Thorsett, S.E., Chakrabarty, D. 1999, ApJ, 512, 288 
Timnies, F.X., Woosley, S.E., Weaver, T.A. 1996, ApJ 
457, 834 

Torra, J., Fernandez, D., Figueras, F. 2000, A&A, 359, 82 
Tsuruta, S. 1979, Phys. Rep., 56, 237 
Tsuruta, S., Teter, M. A., Takatsuka, T. et al. 2002, ApJ, 
571, L143 

Woosley, S.E., Heger, A., Weaver, T.A., 2002, Rev. Mod. 

Phys. 74, 1015 
Yakovlev, D.G., Levenfish, K.P., Shibanov, Yu.A. 1999, 

Phys. Uspekhi, 42, 737 
Yakovlev, D.G., Levenfish, K.P., Potekhin, A.Y. et al. 

2004, A&A, 417, 169 



10 



S. Popov et al.: Population synthesis and neutron star thermal evolution 



Model 


Reference 


TT Cond 


Gaps 


Crust 


T - 
wo BC 


- t 
w BC 


Log N - Log S 


I 


BGV, Fig. 21 
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BGV, Fig. 13 
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BGV, Fig. 15 
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BGV, Fig. 12 
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BGV, Fig. 16 
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BGV, Fig. 14 
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BGV, Fig. 18 
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No 


VIII 


BGV, Fig. 19 
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Yes 


IX 


BGV, Fig. 20 
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Yes 


Yes 
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G05, Fig. 2 
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No 


XI 


G05, Fig. 2 
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B' 
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Table 1. Properties of the selected cooli ng curve s: A - gaps from IXakatsuka fc Tamagakil l)2004f) . neutron gap 
suppressed by 0.1; B - gaps from Yakov lev et alJ l)2004|) . neutron gap suppressed by 0.1; B' - same as for B 
and ^Po proton gap suppressed by 0.5; B" - same as for B, ^Pq proton gap su ppressed by 0.2 and ^Po neutron gap 
suppressed by 0.5; C - - Tin relation fit from BGV; D - - Tin relation bylTsuru3 uM^; E - T^ - Ti„ relation 
from Ya kovlev et alJ 1 20041) and rj — Ax 10"^. The last three entries indicate whether the model complies or not with 
the tests: temperature-age without (wo) or with (w) the additional brightness constraint and LogN-LogS (see text 
and also BGV and G05 for more details). 



